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Possibility for coexistence of a phonon and one-particle spectra in *_ffe at all temperatures in- 
cluding temperatures below the A - point is shown. An approach developed to consider the problem 
of the weakly imperfect Bose gas results in the existence of a gap in the one-particle (one-atom) 
excitation spectrum below the transition point. Thus, the microscopic model of He II turns out to 
differ both from the model initially suggested by Landau and from the Bogolyubov model accepted 
at present. Analysis of the experimental data obtained by measuring enthalpy and the part of he- 
lium atoms condensed in the state with the momentum p = enabled to estimate the size of the 
gap. It turned out to be approximately equal to 0.31 K at zero temperature. 

PACS numbers: 67.00.00, 05.30.-d 



I. INTRODUCTION 

It is generally accepted that the superfluidity (SF) of He II was explained by to L.D. Landau in 1941. Less known is 
the fact that Landau initially assumed the existence of two branches in the excitation spectrum of the system: phonon 
and roton branches. At that, the latter corresponded to quantums of rotational degree of freedom of the liquid and 
was separated from the ground state by an energy gap. 



However, subsequently a viewpoint based, in particular, on the work of Bogolyubov 



ll was accepted. It assumes 



that only one phonon - roton excitation spectrum exists at zero temperature. We think that it is more correct to refer 
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to it as one-particle (one-atom) , because it is believed that the creation and annihilation of phonons cannot lead to a 
change in the number of atoms of the system, while the creation and annihilation operators for the so-called bogolons 
do not commute with the total number of particles in the system. In particular, as a result, the basic Hamiltonian 
of the system obtained by the Bogolyubov transformation does not preserve the total number of atoms. This fact, 
however, is ignored by all the authors: it is considered to be just a consequence of the so-called spontaneous breaking 
of symmetry of the system due to the presence of a second order phase transition. 

This variant of describing the energy spectrum creates a number of questions. On one hand, naturally, no one 
doubts the coexistence of phonon and one-atom spectra at temperatures above the transition point 2]. So, there 
should be a mechanism of disappearance of one of the excitation branches below the transition point. According 
to the conventional point of view the disappearance of the phonon branch occurs immediately after the appearance 
of a macroscopic occupation number of the lower level of the system What experimentalists consider to be a 
spectrum of phonons, strictly speaking, as we mentioned above, is already a one-atom excitation spectrum, although 
it is characterized by a linear relationship between the excitation frequency and the wave vector in the long-wave 
limit. Therefore, a gigantic hysteresis of the transition temperature should exist in the system, because it is obvious 
that the transition point in a boson system with one-atom spectrum corresponding to the experimental phonon - 
roton spectrum of He II should considerably exceed the value corresponding to the parabolic spectrum of the system 
accepted for temperatures above the transition point. In this connection it is pertinent to remind that, according to 
the Landau theory, excitations of the phonon - roton branch form a normal component in the two-fluid model of He 
II. 

In discussing this problem, the systematic research on numerical modelling of the properties of He II in the tem- 
perature range where the system heat capacity considerably differs from the phonon - roton gas heat capacity would 
be extremely interesting and desirable. However, unfortunately, we do not know works in this direction. 

Quite recently the questions, concerning the concept itself of spontaneous breaking of symmetry in the considering 
context, also appeared jsl. The point is that the total number of atoms in the system is not preserved only in case 
of the truncated Hamiltonian, in which only the expressions bilinear with respect to the creation and annihilation 
operators of bogolons are selected. As for the total Hamiltonian of the system, of course, it does not break any 
symmetry. From this point of view the symmetry breaking is rather mathematical in nature and associated with the 
approximations inherent in the method used for calculating the wave functions of the ground state and low-energy 
excitations. But on the physical level, there is no symmetry breaking. If one does not take this into account, a paradox 
Pi results, which until now has not been considered in the works concerning the description of the weakly imperfect 
Bose gas as a qualitative model for the consideration of physical processes in He II in the standard SF theory (of 
course, we assume that the BCS - Bogolyubov theory of superconductivity is obviously associated with the Landau - 
Bogolyubov SF theory). 

Two possible options for further research are possible. On one hand, one can try to prove that solving this paradox 



does not affect tlie final results of calculations of He II within the framework of the microscopic Bogolyubov theory. 
This method seems to be very attractive, because a large number of experimental data were described in the context 
of the standard approach (at least qualitatively). 

On the other hand, we think that it is very interesting to focus efforts on another method: a search of options 
for describing the system, which would maintain strictly the fundamental symmetries of the system at all steps of 
considering it. It is in the search for such an approach that we undertook this study. It will be shown below that it 
suggests a different scheme of the excitations branches in He II. 

II. THE ALTERNATIVE PERTURBATION THEORY IN THE MODEL OF WEAKLY IMPERFECT 

BOSE GAS 



First, note that the result presented in this section was obtained by Bogolyubov |4|1, if to limit to its mathematical 
part. The limitation concerning the mathematical part is not accidental. It is caused by the fact that, when Bo- 
golyubov obtained a gap in the spectrum of quasiparticles, he rejected this solution on the basis of his concept of the 
physical properties of the system. We will discuss this item later. 

So, let us consider spinless bosons (^iJe atoms) with interatomic interaction V (r) as a system located in the box 
Lz ■ We will assume that the boundary conditions of the system are periodical. Then the Hamiltonian of 
the system H in the model of the the weakly imperfect Bose gas {4 [ has the form: 

H ^f + V, 

P,k,q 

where tp — M is the atomic mass of '^He, ft is the system volume, Vq — J V {r) cos (qr) dr , 6+ and bp arc the 
creation and the annihilation operators for "^He atoms in a state with a momentum p. 

Let us calculate the one-atom energy tp of the *He atoms using an equation based on a random phase approximation 
(RPA) H,b^ = tpfop. Let us select first the diagonal part of the Hamiltonian Hdiag , because only this part is 
important for this calculation. For H^iag we easily obtain: 

H,.ag =f+ Y: ^".^^'' Kbib^bp, (2) 

where is the three-dimensional Kronecker symbol. Denominator 1 -t- (5p was introduced in ^ in order to avoid the 
double account of the terms with q = at p = k. 



^diag i '^p 



Taking into account that, obviously, Vq = VLq, equation 
ip: 



tpfep leads to the following expression for 



t,^tp + ^'^Y.-^^^n^, (3) 
k ^ + "p 



where rik is the number of particles in a state with a momentum k. At zero temperature we obtain: 

t +0^^^ (4) 

where c — N/ft is the concentration of atoms in the system. When we obtained ([U, we assumed that at zero 
temperature aU the atoms are in a state with a momentum k = 0, because this approximation takes into account only 
the diagonal terms. 

Obviously, Vb > 0: otherwise, the system would be unstable fl|. On the other hand, is the regular function of 
the parameter q. So, we obtain the gap Aq = Vqc in the spectrum of one-atom excitations of the system. 

When Bogolyubov obtained this result, he thought that it was the result of a wrong description of the system. 
Let us consider in more detail the reason why Bogolyubov considered the appearance of the gap in the spectrum of 
one-atom excitations a wrong solution. 

Analysis of article [4] shows that Bogolyubov expected lim tp/p to coincide with the sound speed in the system and 
thus have a finite value. However, there is no physical principle, from which it would follow that the rate of one-atom 
excitations and the speed of the sound in the boson system should be equal, because the phonons and one-atom 
excitations are two various physical entities. To clarify this point it is sufficient to consider the system above the 
transition point. 

As for the other arguments of Bogolyubov, they are secondary and required to explain why the obtained solution, 
from his point of view, was wrong. These arguments come down to the fact that the perturbation theory becomes 
inapplicable at small excitation momenta. We think that this argument are unjustified, because it is easy to show 
that the appearance of the gap in the one-atom excitation spectrum makes the calculation of the spectrum on the 
basis of the perturbation theory at small excitation momenta more robust. Indeed, the diagonal part of Hamiltonian 
V , which, in fact, is responsible for the appearance of the gap, can be included in the so-called principle Hamiltonian 
of the perturbation theory. Accordingly, if this fact is taken into account, the absolute value of the appearing energy 
denominators turns out to be not less than the value ~ Ap = VpC, which makes the mentioned calculation more 
robust. 

III. PHONONS AND ONE-ATOM EXCITATIONS IN THE PROPOSED APPROACH 

Let us consider an elastic medium, the quantum oscillations of which correspond to creation and annihilation 
operators and /iq. This leads to a new term Hph, to be added to the system Hamiltonian: 

Hph^^sqh+h^, (5) 
q 

where s -is the speed of the sound in the elastic medium thus introduced. Let us name hereafter the quantum 
oscillations of the elastic medium corresponding to the operators and hq seed phonons. 



Let us introduce formally transformations corresponding to the procedure of renormalization of the atom-phonon 
interaction: 



<i^0 q#0 q#0 



2 " V-q i--q rq r q/ p-q- K+q-^k 

k,q#0 

-l^-' E (^r'''''*<+<''>r")^?+q-k/^^'^. (6) 

k^O,q5<iO 

P 

+1^^' E (^r^r'^^' - ^r^vr") hmh+^-^ (7) 

P,k#0 

E (^r'vr - v^r'^vD ^k6+6p_k-q. 

P,k#0 

Operators 6"*" , &, /i^ and h introduced in (l6][7]) are consistent with the standard commutation relations for bosons 
with a precision of ip^. Parameter ip in this approach plays the role of a formal small parameter. 

In RPA approximation we formulate the equation on the transformation parameters such a way that the system 
Hamiltonian in terms of operators , b, and h would not include terms of the structure h^h^h and hb^b. This 
equation will be of the form: 

- V^Jg + (tp_q -tp+ Sq) + Vq \Jg\' + jl^Pq'l* = 0, (8) 

where 

J, = f7-i^(^P+'i-^P)np, (9) 
p 

and Up is the thermodynamic average of the number of atoms in the state with a momentum p. In ([8]) the numerical 
coefficient with terms ^ Vq\Jq\^, ~ VqJq was chosen in such a way that the spectrum calculated in the procedure 
of renormalization (see (|15p ) coincides with the spectrum of sound obtained on the basis of analysis of the spectral 
dependence of the system dielectric function [5] . Actually, in this case, we simply use the opportunity of choosing the 
order in which the terms are considered in a situation where the system is missing a real small parameter. 

A nontrivial solution of equations (I5][7]) under the assumption that parameter ip^ is real- valued is given below: 

^ JMtp-tp^q + sq)^ ^^^^ 



At that, ujq is consistent with the foUowing equation: 

=0. (12) 

Thus, the dependence of ujq gives a sound spectrum obtained by analyzing the dielectric function of the system. 

Kinetic terms /fp/" contributing to the rcnormalization of energy of the seed phonons calculated with the use of 
RPA are of the form: 



= E 'lK''^ + Hv + Hu (13) 
q 



where 



q q q 

Ht = J2 ((Wq -tp~ si) k^'^l' + (ip-q -tp+ kq|') "p^^q^' 



(14) 



p,q 

+ J2 {tp - ip + q) y^rV^q^P^q '^iq 

P.q 

+ E (ip - tp + q) ^P+^^%nphqh_q. 

p.q 

For clarity we will first take into account only the contribution of Hv to H!^™ , which, actually, as is shown below. 



corresponds to the calculation of the spectrum in the long- wave limit. Diagonalization of the operator -ff^J" in this 



ph 

-1 1 1 r7 o ^- 1 /-111 /-if /-iT-idyo ^-/~lT■ P, 

case using the Bogolyubov transformation leads to the following expression for the rcnormalized energy of the seed 
phonons Sq. 

Eq = Wq, (15) 

As is expected, the spectrum of the seed phonons coincides with the expression obtained for the spectrum of sound 
vibrations by calculating the dielectric function. In particular, for the Bose gas at zero ternperature coincides 
with the well-known expression obtained by Bogolyubov for the spectrum of quasiparticles l| with neglect of the 
difference between the squared absolute value of the so-called Bogolyubov c - number and the full number of atoms 
in the system, i.e. we have: 



tq(tq + 2Aq). (16) 

The resulting expressions for the rcnormalized quantities in the suggested approach should be taken in the limit 
s 0. Indeed, in this case Hph does not add terms to the system Hamiltonian. So, the system of seed phonons only 
gives a set of states, which are used for the quantization of its oscillations. In particular, for the additional term of 
interatomic interaction SV following from the terms of Hamiltonian T and Hph, we can obtain: 

sv ^ ^n-' ^y^'^'bUAJ'^K, (17) 



where 



1 '^q 



(18) 



Considering equation (jlSp leads (in the limit p,k 0) to the conclusion that rcnormalization docs not contribute to 
the size of the gap in the one-atom excitation spectrum. 

Now let us take into account also the contribution of Ht in the spectrum of the phonons for the case of zero 
temperature. A simple calculation gives the following expression (compare with (jl5p ): 



In fact, one can say that in the system one more expansion parameter appears, which is equal to iq/2Aq. It can be 
assumed, therefore, that the spectrum of phonons ends at the values of wave vector q corresponding to the condition 
tq ^ 2Aq. Recall in this context that the spectrum of phonons is also limited at the shortwave region when the 
Boltzmann distribution for the occupation numbers of the atomic energy levels is determined. Besides, at high 
temperatures T > Aq/4 equation l|12p has no solution for cjq. 

The suggested method makes it possible in principle to obtain the properties of the system with any accuracy, if one 
does not take into account the arbitrariness in the order of consideration of the terms appearing in the rcnormalization 
procedure. This arbitrariness is caused by the lack of a real small parameter in the system. We plan to eliminate it 
by comparison with the results of alternative methods of system's investigation. 

Of course, not all the renormalized values will be calculated so simply as the correction to interatomic interaction 
energy SV^ because transformation parameter ip^ has a singularity at s — > 0. At the same time, it seems that this 
problem is not of a fundamental nature. It is sufRcient to increase the accuracy of transformations ([6][7]). At that, the 
appearing singular expressions will be reduced leaving only the finite values in the calculation results. 

In summary, one can assert that at this stage of development of the suggested approach no reasons can be seen why 
the coexistence of the two branches of the excitation spectrum of the system may be doubted at temperatures below 
the transition point. It is quite another matter that the results obtained in the context of the suggested approach 
should be confirmed by other methods. Of course, in this case we mean a description of the system with the use of 

n . 

the technique of Green functions modified as compared with 6] . This technique should not break the fundamental 
symmetries of the system. We plan to perform such consideration in near future. 

IV. BOSE GAS THERMODYNAMICS OF ONE-ATOM EXCITATIONS WITH INTERACTION 

Now let us consider the case of finite temperatures assuming that they are sufficiently low to carry out the condition 

^ Vk-prik « VpiVe, (20) 



8 

where is the full number of atoms in the excited states. In this case the simplification of ([3]) gives the following 
expressions: 

?o = 2Ao - kAo, 

(21) 

ip = tp + 2Ao + (Ap - Ao) {at p 7^ 0} , 
where k, is the fraction of atoms condensed in the ground orbital. 

Let us assume that in the case of thermally activated atoms it is possible to use the parabolic dispersion law with 
a mass Afe, that is 

tp«2Ao + ^ (22) 

at p 7^ 0. This simplification is possible, because Vp — Vq ^ at p ^ 0. 
As for the system energy E , using the expression for Hdiag gives: 



E = Y, <P"P + ^ E "P ("P - 1) + \^'^ E + ^k-p) ripn^ (23) 

P P P5^k 

Again, assuming that the temperatures are sufficiently low to satisfy the conditions 



E V^,_pnpnk~Vb E npTi^^^VoNl, (24) 

k^O.PT^O.pT^k k5<iO,p#0,p#k; 



it is possible to calculate energy e per one atom of the large system: 

,2 
2 



e = Ao ( 1 - ^ ) + (1 - «;) (t), , (25) 



where {t) ^ is the average kinetic energy of the excited atoms of the system: 



(*)e=^e"'E^p(^p)' (26) 
Pt^O 



For average occupation numbers (rip) the Bose-Einstein distribution function is applicable. Thus, we have: 

, , 1 



(27) 



^ exp ((p2/2Me + kAo) /T) - 1 ' 

because we assume, as usual, that the value of the chemical potential of the system is very close to io- 
Pass from sums to integrals it is possible to obtain an equation for the fraction of condensed atoms: 

47r2c J exp {kAq/T + x) - 1 ^ ^ 

and an expression for (i)^: 

r x^/xdx 
m ^e^ J cxp(KAo/T+a;)-T 

We - ^rf, - ■ (29) 

J gxp(kAo/T+x)-1 

Accordingly, it becomes possible to compare experimental data with calculations based on the proposed scheme and 
to estimate the gap Aq and mass A/e of the excited atoms of the system. 



9 



TABLE I: Basic data and calculation results for a temperature range 1.5 — 2.1K 





1.5 


1.6 


1.7 


1.8 


1.9 


2.0 


2.1 


Utot, cal/g [7] 


0.071 


0.112 


0.176 


0.252 


0.351 


0.484 


0.662 


Uph, cal/g [8] 


0.046 


0.068 


0.095 


0.129 


0.170 


0.219 


0.276 


K, % [9] 


94.33 


89.83 


83.92 


76.12 


66.90 


55.79 


41.84 


Uat, cal/g 


0.025 


0.044 


0.081 


0.123 


0.181 


0.265 


0.386 


A,K 


0.19 


0.08 


0.04 


0.01 
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0.43 


0.52 


0.63 


0.70 









V. DETERMINATION OF MODEL PARAMETERS BASED ON EXPERIMENTAL DATA WITH THE 
USE OF ENTHALPY AND THE FRACTION OF THE SF COMPONENT IN HE II 



The calculation idea consists in the possibihty of using two equations (1251 128|) in order to determine two parameters 
Aq and 7 = M^/M at known values of temperature T, energy e, and fraction of atoms condensed at the ground orbital 



K. 

The value of energy e per one atom of helium can be obtained from the data of Kapitsa on the temperature 
dependence of enthalpy [vj, because under the conditions of an isobaric experiment it is very close in value to the 
temperature dependence of the internal energy of the system. Indeed, pressure in the experiments of Kapitsa does 
not exceed 25 atm. So, the difference between the enthalpy of helium at zero and finite temperatures is very close 
in value to the difference of the corresponding internal energy values, because density changes in He II caused by 
temperature are very small. 

According to the proposed concept the total enthalpy of the system determined in the experiment is a sum of 
two contributions: phonon - roton and one-atom excitations. Now we have a task of extraction the contribution of 
one-atom excitations from the total energy of the system. 

The phonon contribution to the internal energy of the system is connected with the temperature dependence of the 
phonon heat capacity by the relation: 

T 

UphiT)^ JCphit)dt. (30) 





At that, there is the well-known expression 8] for the heat capacity Cph {t) : 

27r2, /fcst^^ 



pWfJ,rkBt/2 ( Al Ar 3\ f Ar' 



where ks - is Boltzmann constant. Nowadays experimental values of the phonon spectrum parameters used in ([31 
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T, K 



FIG. 1: Temperature dependence of one-atom excitation energy Se calculated from atomic energy Uat per one helium atom 
(upper curve). Besides, a calculated curve for the dependence 5e{T) on the basis of (|34p at selected values of Aq and 7 
determined by (|36|l is presented (lower curve). All data are given in Kelvins. 



(32) 



are as follows: 

S = 2.38 • 10"* cm/sec, 

= 8.6 K, 
po = 1-8 • 10^ n/cm, 
fir ^ 0-13 M. 

Obviously, within the limits of a simple model of two branches of excitations one-atom excitation energy Uat (T) can 
be obtained by simple subtraction of the phonon-roton contribution from the total energy of the system Utot (T): 



Uat (T) = Utot (T) - Uph (T) . 



(33) 



Temperature dependence of the fraction of condensed atoms in the ground state k (T) can be directly obtained from 
the results of an experiment by Andronikashvili (measuring the period of axial-torsion oscillations of a vessel filled 



with He II [9], or that by Peshkov (measuring the speed of second sound [10|). Both experiments give very similar 
dependencies n{T). It is relevant to note that the fundamentals of the phenomenological theory of He II, including 
the possibility of a two-fluid description of the system, cannot be doubted. 

Table 1 presents the basic data for calculation of Uat (T) , as well as the results of calculations of other parameters 
that will be determined below. Fig. 1, respectively, presents the dependence of Se (T), value equal to energy Uat (T) 
calculated per one atom of helium. Note here that the energy Se, obviously, should be equal to energy e, counted 
from its own value at T = i^T, i.e. (compare with ((25|) ) 

(34) 



Se = Ao (1 - K^) /2 + (1 - k) -fT 



f : 



cxp(KAo/T+a:)-l 



r ^Jxdx 
J cxp(KAQ/T+a:)-l 



Dependencies A (T) and 7 (T) presented in the table are determined with the use of a self-consistent system of two 
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0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



80 



85 



90 



95 



100 



k, % 



FIG. 2: Calculated dependencies A (k) (solid curve, A are given in Kelvins) and 7 (k) (dashed curve) 



equations (compare with 



6e (T) = A (T) (1 - (T)) /2 



I 



'^^^) = ' A^c J 



cxp(k(T)A(T)/T+x)-1 

/\fxdx ' 
oxp(k(T)A(T)/T+i)-1 

' \Jxdx 



(35) 



exp K (T) A (T) /T + a; - 1 



The obtained dependencies A (T) and 7 (T) are the basis for obtaining the parameters Aq and 7 of the model. 
Indeed, when we deduced formula ([25|) we made the assumptions (120 1 124 |) ). which are realized the better the lower 
is the system temperature. On the other hand, it is difficult for us to estimate the temperature starting at which 
assumptions (j201l24|) are realized with acceptable accuracy, because there are no reliable data on interatomic potentials. 
In this situation in order to obtain parameters Aq and 7 it is necessary to extrapolate dependencies A (T) and 7 (T) 
to the zero temperature range, where these assumptions are realized perfectly. 

On the other hand, using the dependency k (T) enables to obtain the dependencies A (k) and 7 (k) shown in Fig. 2. 
At that, it becomes possible to extrapolate the obtained functions to the range of k = 1. This also corresponds to the 
zero temperature case. This option seems to be more attractive, because the extrapolation becomes non-dimensional, 
and its value of « 5% seems to be small. Extrapolation parameters Aq and 7 obtained in this case have the following 
values: 

Ao = 0.31 K, 

(36) 

7 = 0.28. 

It is obvious that dependency 8e (T) , calculated on the basis of ([M)) at the values of parameters Aq and 7 selected 
in ([55)) should approach dependency Se (T) calculated on the basis of ([55)1 at decreasing temperature, because in this 
case correlations (1^ [M| are realized better. Fig. 1 visualizes this statement. 
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VI. CALCULATION OF THE DIELECTRIC FUNCTION OF THE SYSTEM 

When calculating the system response to an external action one should keep in mind the requirement of gradient 
invariance of the model. The below calculation takes into account this circumstance by using the well-known fact 
used for calculating the Hartree dielectric function: if Hamiltonian T is taken as the basic Hamiltonian of the system, 
it is possible to calculate the response on the external potential with use of the standard theory of excitations. The 
Hamiltonian of this external potential tj has the form 

U = j U{v,t)h+{v)b{v), (37) 

where fo"*" (r) and 6 (r) are creation an annihilation operators in coordinate representation. In this case the gradient 
symmetry of the system is not broken. 

As one goes to a new basis of wave functions, in order to solve this problem in the context of the perturbation 
theory, a number of terms appear that have the same structure in RPA as the right hand side of equation (j37p . 
This circumstance is taken into account by introducing effective potential Ueff (r,t). At that, for Fourier transforms 
U (r,t) and Ueff (r,t) the following correlation exists: 



where Eh (q, is the Hartree dielectric function: 



The presented procedure does not take into account the residual Hamiltonian Hres, for which we have: 

Hres--V + U-Ueff, (40) 

where 

Ueff = j Ueff {v,t)h+{v)h{v)dv. (41) 

In the new basis of wave functions, obtained after solving the Hartree problem of the perturbation theory, the matrix 
elements Hres are invariant under the gradient transformation. Further let us restrict consideration to a simple case 
of the diagonal terms. 

This consideration leads to a shift of the system energy levels. It turns to be important for determination of the 
occupation numbers rik and nk+q in (|39p : calculating them one should use the Bosc-Einstcin distribution function, 
in which one-atom energies ik and tk+q are used as energy levels. 

Note one more feature of the Hamiltonian Hres- Obviously, the response of the system on the Hamiltonian Ueff 
described by the dielectric function at zero temperature corresponds to the response of an ideal SF fluid. At the same 
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time, non-diagonal, time-dependent terms of the Hamiltonian Hres can lead to dissipation of energy due to transitions 
from the ground state to the excited ones if the frequency of the external action exceeds the gap energy in the one- 
atom excitation spectrum. Here it is reasonable to make a most general remark: the appearance of time-dependent 
non-diagonal terms of the Hamiltonian in the procedure of calculating the system response to an external action for 
the models using gapless one-atom excitation spectra can mean the appearance of an energy dissipation channel in 
the system, even if it satisfies the SF principles from a formal point of view. 

On the other hand, the hypothetical possibility of the existence of a gap in the one-atom excitation spectrum of 



SF state has been suggested before 
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12l |. For such models it is possible to observe directly the gap size values. For 
example, if voltage is applied to a metal needle, then near its tip a force appears that exerts influence upon the helium 
atoms. Therefore, the conditions of the experiment proposed earlier in very general terms 1J|, become better defined, 
because now we know the frequency interval, in which a thorough investigation of the phase and amplitude-frequency 
characteristics of the chain which includes a sensor in the form of the above needle should be carried out. 



VII. CONCLUSION 



This paper suggests the SF model of He II preserving the fundamental symmetries of the system at all temperatures 
including those below the transition point. The structure of energy levels appearing in the model is determined by 
the possibility of coexistence of phonon and one- atom spectra below the transition point. 

Numerical simulation showed satisfactory agreement of the model with the experimental data below the transition 
point. In particular, we succeeded in estimating the value of the gap in the one-atom excitation spectrum, which 
turned out to be approximately equal to 0.31 K at zero temperature. We suggest a method of direct observation 
of the gap in the one-atom excitation spectrum on basis of the data on energy absorption in the system under the 
influence of an external oscillating force. 
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